Two-dimensional imaging of gauge fields in optical lattices 
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We propose a scheme to generate an arbitrary Abelian vector potential for atoms trapped in a two- 
dimensional optical lattice. By making the optical lattice potential dependent on the atomic state, 
we transform the problem into that of a two-dimensional imaging. It is shown that an arbitrarily fine 
pattern of the gauge field in the lattice can be realized without need of diffraction-limited imaging. 



Recently, many-body systems of trapped atoms have 
offered a new avenue towards understanding strongly cor- 
related matters In this context, trapped atoms in 
gauge fields have attracted much attention [2|-[5|. In a 
two-dimensional (2D) confinement, such systems would 
exhibit intriguing phenomena, such as the fractional 
quantum Hall effect and the anyonic statistics, which 
have opened up a new era of modern condensed matter 
physics 0,0. 

The aim of this paper is to propose a concrete frame- 
work to realize the following Hamiltonian in an optical 
lattice system 



3,k 



C j,k+l C 3, k ' 



H.c. 



(1) 



which describes a tight-binding particle in the lowest 
band of a 2D square lattice in a spatially slow-varying 
magnetic field B(r) = V x A(r) 0|. Here, the positions 
of lattice sites are = ro(jx + ky) with ro being the 
lattice spacing, c 3 ^ denotes a particle annihilation oper- 
ator at the site, and = (27r/$ ) J*** 1 '* A • dl with 
$o = h/e being the magnetic flux quantum. We have 
chosen a gauge such that the vector potential is written 
as A(r) = A(x, y)x. While the lattice structure gives 
rise to intriguing physics in its own right, such as Hofs- 
tadter butterfly [9[ , this model Hamiltonian also restores 
the fractional quantum Hall physics in the continuum 
limit when a strong on-site interaction is considered in 
addition [5|. Such an optical lattice system would thus 
provide an ideal testbed to study the physics in the pres- 
ence of gauge fields, along with established techniques for 
its microscopic control and measurement [Io|, [TlJ • 

Before proceeding further, we highlight the underly- 
ing idea and important features of the present work. 
Differently from earlier schemes in this context 0, Q, 
our scheme is devised to make it feasible to realize an 
arbitrary A(r), which enables the generation of quasi- 
particles/quasiholes by adiabatically changing A(r) 
This may open up an intriguing possibility of directly 
and microscopically inspecting the fractional charge and 
statistics of the quasiparticles, as well as other modes of 
measurement discussed in earlier literature [3[. Further- 
more, we achieve this goal along with overcoming the 
earlier experimental problems, such as the stability issue 
or complicated experimental setup d, Q . Unlike most of 



the earlier schemes based on engineering Raman-induced 
hopping, we simply make use of the ordinary tunneling in 
an optical lattice. We employ the state-dependent lattice 
potential that naturally comes about by adjusting the de- 
tuning and polarization of the trapping laser [12| . On top 
of this simple setup, we show that the spatial phase pro- 
file of a static driving field, which induces on-site Raman 
transition of atoms with different phases, is mapped into 
the Abelian gauge field A(r). Our scheme should thus 
be well within state-of-the-art technology and compati- 
ble with currently prevalent setups based on alkali-metal 
atoms. Moreover, as we use the atomic dark states de- 
coupled from the Raman field, the recoil heating is min- 
imized. While the diffraction limit seems problematic in 
imaging an arbitrary phase pattern site-by-site on the 
lattice, we show that a 2D array of only moderately well 
focused beams, each centered on a lattice site, in fact 
suffices for any A(r), aside from the available technology 
of subdiffraction imaging [13[. 

We first present our (abstract) Hubbard model, 
putting off the detailed explanation of its realization un- 
til later. Let us consider a situation where two hyperfine 
levels | a) and \b) of an atom experience different opti- 
cal lattice potentials, while the two lattices coincide. In 
particular, let us assume the atoms in \a) (|6)) hop only 
in the x (y) direction. The Hamiltonian for this state- 
dependent hopping can be written as 



Hi = -J J2( a ]+hk a j,k + &3, fc+ i&i,fc + H.c), 

3,k 



(2) 



where a^k (pj,k) denotes the annihilation operator for 
atoms in \a) (\b)) and for simplicity the hopping rate J is 
assumed to be the same for each direction. On top of this, 
Raman fields for transition \a) \b) are applied with 
position-dependent phases fy^. Note that these phases 
are determined unambiguously by the difference between 
the phases of two Raman fields. The Hamiltonian for the 
Raman transitions is written as 



'J2( a lk b j 

3,k 



-i<Pj, k + h.c.) 



(3) 



where uj is the (real-valued) Raman transition rate as- 
sumed to be the same for every lattice site. For con- 
venience, let us perform a local gauge transformation 
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FIG. 1. The smallest loop in a square lattice (a) and a bilayer 
lattice we consider (b) . In (b) , the upper and lower layers trap 
atoms in \a) and respectively. Both loops enclose the same 
number of magnetic flux quanta (mod 1) (6j,k — #j,/c+i)/27r. 



dj y k —> aj y ke z ^'' fc . The single-particle Hamiltonian H s = 
Hi + H r then reads 

H * = - J E( fl kfc fl ^ e " i,fe + b lk+ihk + h.c.) 

3,k 



(4) 



where Oj^ = <fij+i,k — </>j,k- Note that 0^ can be chosen 
arbitrarily and independently for each pair of sites. We 
are going to choose Oj^ as defined in Eq. (pQ) (conversely, 
the laser phases determine the vector potential A(r) of 
the simulated system) . Note that the smallest loop of this 
system as shown in FIG. H^b), which is not 2D, encloses 
the same amount of magnetic fluxes as that in FIG.[TJa). 

A further progress is made by performing a trans- 
formation such that Cj^ = ~^( a j,k ~ and dj : k = 
~^( a j,k + fy/,fc)- Eq. (HJ then reads H s = Hq where 



H.= 



J oE( c j+i/i^ e%fc + c lk+i c 3\k + H.c.) 
- JoY,( d ]+hk d j,ke i<>j ' k + 4,k+i d i,k + H.c.) 
+ ujJ2 { d ),k d i,k - c) tk Cj, k ) , (5) 

j,k 

Hi=- JoY,( c l+i,k d 3 ,ke ie - k + d) +ltk c hk e^ + H.c.) 

j,k 

+ J ° J2( c ik+i d 3,k + dt k+1 c jik + H.c), (6) 

3\k 

where Jo = J/2 is chosen in accordance with Eq. (pQ). 
It can be seen that the energy spectrum of Ho alone is 
divided into two bands centered at ±gj, each corresponds 
exactly to that of our desired Hamiltonian (pp), whereas 
Hi only contains terms that give rise to transitions be- 
tween the two bands. These transitions are energetically 
costly for 2uj ^> Jo, in which case they are suppressed 
and hence Hi can be treated as a perturbation. This be- 
comes clearer when we consider a uniform magnetic field 
Oj^k = 2nak, where a is the number of magnetic flux 
quanta passing through the lattice cell. In FIG.[2j we plot 
the energy spectrum of Hamiltonian (jl]). The horizontal 
axis represents the energy and the vertical axis a, where 



rational numbers a = p/q with p, q coprime integers and 
q < 50 are taken. As the width of each energy band is 
4 J, two energy bands are separated when uo > 2 J. As 
is expected, for a large uo (e.g., see FIG.[2fd)), the origi- 
nal shape of the Hofstadter butterfly is restored in each 
band, albeit slightly deformed due to the perturbation. 
This deformation gets smaller for larger uj. 

For a system of N atoms, one can choose the lowest 
energy band centered at —Nuu, wherein all the particles 
are approximately in Cj^ modes. It is noteworthy that 
Cj^ particles are in dark states when the corresponding 
two Raman fields have the same Rabi frequency. As these 
particles are decoupled from the Raman fields, they are, 
in principle, not influenced by the spontaneous emission. 
In order to see that this system can produce the fractional 
quantum Hall physics, we again employ a uniform mag- 
netic field Oj^k = ^irak and consider the following on-site 
interaction Hamiltonian on top of the above single-atom 
Hamiltonian: 

3,k 

where for simplicity all the interaction rates are assumed 
to be the same as U . As an example, we have numerically 
diagonalized the Hamiltonian for two atoms in a 8 x 8 lat- 
tice, where a periodic boundary condition with respect to 
the magnetic translation is taken so as to account for a 
large lattice system [14|. We have taken U = uj = 10 J 
and a = 1/16, hence the filling factor is v = 1/2. Due 
to the torus geometry, this system has twofold degen- 
erate ground states. We have obtained the two lowest 
energy eigenstates p M (/i = 1, 2) by tracing out the inter- 
nal state. The purities of these states are Tr(p^) = 0.994 
and Tr(p^) = 0.993, which means the internal degree of 
freedom is almost decoupled from the motional degree 
of freedom, and the number of particles in Cj^ modes is 
c ],k c 3,kP/j,) ~ 2.000. These are in good agreement 
with the discussion above. We also define a projector to 
the subspace spanned by the two Laughlin ground states 
on the torus P L = (*f | + |*2 ) (*2 | and calculate 
the overlap of the above numerical ground states to this 
subspace, which is found out to be Ty(PlPhPl) = 0.999. 
We also consider second-nearest-neighbor hopping in case 
the lattice potential in the direction of hopping should be 
set relatively shallow. As this does not destroy the struc- 
ture associated with a, the ground state is not expected 
to be altered significantly. Our numerical calculation in- 
dicates that for the second-nearest-neighbor hopping rate 
set to be J/ 10, the overlap of the ground state is still as 
high as 0.998. 

The remaining question is how to realize Hamiltoni- 
ans (|2]) and (|3]). To be concrete, we introduce below a 
particular setup based on alkali- metal atoms, generaliz- 
ing the conventionally used one-dimensional (ID) setups 
[12j . However, the idea would be applicable to different 
atomic species or laser configurations. 
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FIG. 2. The energy spectrum of the system for different parameter choices. In (d), only one of the two bands, with the energy 
ranging from —12 J to —8 J, is plotted, which resembles the original Hofstadter butterfly shown in (e). 
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FIG. 3. (a) The fine structure of alkali-metal atoms, where 
the transitions by <j_ polarization, equivalent to those by <j+ 
with mj — >> -mj, are omitted for brevity, and (b) a laser 
configuration for the lattice potential in x direction. 



Hamiltonian (j2]) requires the lattice potential to de- 
pend on the atomic state. This can be achieved by 
exploiting the polarization dependence of atomic dipole 
transitions. Let us first consider the hopping in x di- 
rection. The basic idea is as follows [12|. Consider the 
fine structure of alkali-metal atoms and suppose that the 
frequency of the trapping laser is chosen between those 
of Di (S1/2 P1/2) and D2 (S1/2 P3/2) transitions, 
as shown in FIG. [3ja). If the trapping laser is in <j+ 
polarization, atoms in |+) would experience a negative 
ac Stark shift V+ as the laser is red-detuned from the 
D2 transition. On the other hand, atoms in |— ) would 
experience both negative and positive ac Stark shifts 
due to the D2 and Di transitions, respectively. When 
these two Stark shifts are summed, the resulting en- 
ergy shift V- can range from negative to positive val- 
ues, depending on the laser frequency, while V+ is still 
negative in any case. The ratio V+/V- can thus be 
adjusted to a great extent. The energy shift for each 
hyperfine level is now obtained as a linear combination 
of V±, determined by the Clebsch-Gordan coefficients. 
As an example, let us take \a) = \F = = +1) 

and I b) = \F = 2, mp = +1) among the ground hyper- 
fine levels of 87 Rb. The ratio between the correspond- 
ing lattice potentials is then given by Vb(x)/V a (x) = 
(3V + + V-)/(V+ + 3 VI). For example, suppose the fre- 
quency of the laser is chosen such that V + = — 7V_, 
leading to Vb(x)/V a (x) = 5 (note that this can be done 
when the frequency is red-detuned compared to that cor- 



responding to V- = 0). In ID cases, the hopping rate 
is proportional to E r (Vo/E r ) 3 / 4 exp(— 2^/Vq/ E r ), where 
E r = h 2 k 2 /2m is the recoil energy and Vo is the poten- 
tial depth [lOj. If we take moderately V a (x) = —bE r 
at the minima, the above choice leads to a hopping rate 
for I b) being 0.013 times that for \a) in x direction, well 
approximating our model. 

In an analogous fashion, one could use <r_ polariza- 
tion in y direction. In a 2D geometry, however, this 
cannot be done without introducing tt polarization when 
the beam in x direction is chosen to be <r + polarized. 
When the rc polarization is involved, as well as the ac 
Stark shift, a Raman transition can also take place, as 
can be seen in FIG. [3fa). Although the Raman transi- 
tion could be suppressed by applying a magnetic field, 
it would be advantageous to exploit the large fine split- 
ting to obtain a stronger potential, as we discuss below. 
Let us take the convention that a± = -^(x ± y) and 
tt = z. In FIG. [3fb), we depict the laser configuration 
in x direction. The lattice potential in y direction can 
be created analogously. In x direction, two tilted stand- 
ing waves with wavevectors ±fci = ±k(cosrjx — sin r\z) 
and ±&2 = ±k(cosr]x + sin^i) and polarizations e\ oc 
V^sin?]^ + cos r]z and £2 oc \/2 sin 770+ — cost?£, re- 
spectively, are applied (note k\ • e\ — k,2 • £2 — 0). 
Choosing the phases appropriately, the electric field 
can be written as E(r) = a+E+(r) + irE^fr) with 
E+(r) ex V2 sin r] cos(kx cost]) cos(kz sin r]) and ^(r) oc 
cos 77 sin (/cx cost?) sin(kz sinr]). Note that E+(r) and 
E n (r) are ir/2 out of phase in both directions. The min- 
ima of the potential for <r+ (oc — \E+(r)\ 2 ) correspond to 
the maxima of the potential for tt (oc — ^^(r)] 2 ), and vice 
versa. If we choose a 2D plane formed, e.g., at z = with 
a strong confining potential in z direction, the influence 
of the tt polarization to the lattice potential is negligible. 
Furthermore, denoting by W±(r) the (real and symmet- 
ric) Wannier functions for state the Raman transi- 
tion rate oc \J <ir,E + (r) J E 7r (r)W + (r - Yj^)W-(r - Tj^)\ 
always vanishes for every site because the electric 
field part is an odd function around r^. Note that the 
lattice spacing is 1/ cos 77 times larger than usual A/2 be- 
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cause the lasers are tilted. For the same Vb, this decreases 
the hopping rate, hence the energy scale of the system. 
Whereas this can be mitigated by decreasing 77, smaller 
r] results in a shallower potential. This trade-off relation 
would be important in choosing 77 in experiments. 

Finally, we discuss the realization of Hamiltonian (j3j). 
For the hyperfme levels \a) and \b) chosen above, the 
Raman transition in Hamiltonian (|3]) can be achieved 
by applying two <r+ polarized Raman fields with differ- 
ent frequencies. As the 2D lattice is on the x-y plane, 
this can be naturally done by applying the light in z 
direction. Suppose the spatial profiles of the Raman 
fields are given by Q{ a ^ b y(x, y)e lkz . The Raman tran- 
sition rate at site is then given by ujj^ l ^ j,k oc 
/ Ql(x,y)Q b (x,y)W a (x - jr ,y - kr )W b (x - jr ,y - 
kr$)dxdy, where W a (x,y) and W b (x,y) are the Wan- 
nier functions for \a) and |6), respectively, and the 
z dependency is assumed to be decoupled. Suppose 
n* a (x,y)Q b (x,y) ~ IOqI 2 ^'^. We require <j){x,y) to 
be slowly varying compared to W{ a ^y(x,y) so that ujj^ 
to be the same for every site. <pj^ is then approxi- 
mately given by (f)(jro, kro). Note that any nonlinear 
term in <j)(x, y) can lead to non- vanishing magnetic flux 
(Oj,k — 0j,fc+i) 7^ 0- For example, computer generated 
holography would be relevant for adjusting <j>(x , y) arbi- 
trarily as far as its resolution limit permits |15| . 

A truly arbitrary magnetic field would be possible 
when 6j y k (or (j)j,k) can be controlled site by site. At 
first sight, this seems to require diffraction-limited imag- 
ing. For example, this is the case for a uniform magnetic 
field (j)j y k = 2itajk (mod 2ir) in the sense that for lat- 
tice sites with j = l/2a, the phase = (—l) k 7r should 
be flipped at every adjacent site. As will be shown be- 
low, however, this seemingly necessary requirement of 
diffraction-limited imaging is in fact not essential in our 
scheme, although it should also be noted that subdiffr ac- 
tion imaging technologies are already available [13]. For 
convenience, let us represent the sites (j, k) with a sin- 
gle index £ or A. We apply at each site £ a moder- 
ately well focused Raman field for \b) with spatial pro- 
file (j0^e l ^^A(r — r^), where A(r) is a common normalized 
mode function, e.g., a Gaussian function, while the other 
Raman field for \a) is applied globally and uniformly. The 
Raman transition rate at site A is then given by the sum- 
mation of all the contributions uje lc ^ x oc J^T^c^e*^, 
where T\£ = J Air - r^)W a (r - r x )W b (r - r x )dxdy. For 
well focused beams, T\% would be non- vanishing only for 
small \r\ — r^|. The problem now reduces to determining 
a proper set of uj^ and (j)^ to obtain the desired field pro- 
file corresponding to the given set of It is easily seen 
that this can be done by calculating the inverse matrix 
of T for any vector uje l ^ x , i.e., any A(r). A similar idea 



has been used in a different context in Ref. [l6| . 

In summary, we have shown that in a state-dependent 
optical lattice, a light field itself can play the role of a 
gauge field. This is an intriguing situation where the 
phase gradient of light solely produces a torque to atoms. 
While this has been known in a different context [I?}, it 
will be an interesting future work to clarify their analo- 
gies. 
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